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Abstract. Let be a Lipschitz domain in R™ n > 2, and L = divAV be a second 
order elliptic operator in divergence form. We establish solvability of the Dirichlet 
regularity problem with boundary data in H ,p (9f2) and of the Neumann problem 
with L p (dil) data for the operator L on Lipschitz domains with small Lipschitz con- 
stant. We allow the coefficients of the operator L to be rough obeying a certain 
Carleson condition with small norm. These results complete the results of [S] where 
L p (dil) Dirichlet problem was considered under the same assumptions and [5] where 
the regularity and Neumann problems were considered on two dimensional domains. 



1. Introduction 

This paper continues the study, began in [5], of boundary value problem for second 
order divergence form elliptic operators, when the coefficients satisfy a certain natural, 
minimal smoothness condition. Specifically, we consider operators L = div(v4V) such 
that A(X) = (a,ij(X)) is strongly elliptic in the sense that there exists a positive 
constant A such that 

for all X and all £ G W 1 . We do not assume symmetry of the matrix A. There are a 
variety of reasons for studying the non-symmetric situation. These include the connec- 
tions with non-divergence form equations, and the broader issue of obtaining estimates 
on elliptic measure in the absence of special L 2 identities which relate tangential and 
normal derivatives. 

In [H], the study of nonsymmetric divergence form operators with bounded mea- 
surable coefficients was initiated. In [11], the methods of [14] were used to prove 
results for the elliptic measure of operators satisfying (a variant of) the Carleson mea- 
sure condition. 

Our main result is that under the assumption that 
(1-1) ^(Xy 1 (osc B{x ,8(x)/2)aij) 2 

is the density of Carleson measure with small Carleson norm then the Dirichlet problem 
for operator L with boundary data in H 1,p (dVL) is solvable. We also obtain the same 
result for the Neumann boundary value problem with L p (dfl) data. 

Let us recall that the paper [5] considered the L p (dQ) Dirichlet problem under the 
same assumptions on the coefficients. It turns out that the regularity and Neumann 
boundary value problems are considerable more difficult and the progress on these 
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problems under assumption (11.11) has been fairly slow. Only recently these two prob- 
lems have been resolved on two dimensional domains in [5J. However, the proof relies 
fundamentally on the fact that the domain is two dimensional. In particular we are not 
currently able to generalise the concept of a conjugate solution to higher dimensional 
cases in our context. We should note, however, that such generalisations have been 
carried out in other contexts (see Theorem 9.3 of [1] and Section 1 of [16] which puts 
the result in context). Our paper relies on several recent advances, in particular [3] 
and [1] where a better understanding of the Dirichlet and regularity boundary values 
problems was obtained including the "duality" between the solvability of the Dirich- 
let boundary value problem and the regularity problem for the adjoint operator. In 
particular, this allows us to focus on solvability of the regularity problem for p = 2 
and obtain solvability for values of p with no additional effort. We are able to avoid 
the use of the p-adapted square function introduced in [3] and which was essential in 
establishing solvability of the L p Dirichlet problem. The solvability of the Neumann 
problem is trickier as there is no appropriate analogue of the results in [3]. We over- 
come this by using the solvability of the regularity problem, which we established here, 
and induction (see Lemma [6. 2p . 

Operators whose coefficients satisfy small or vanishing Carleson condition (II. ip arise 
in several contexts. For example, consider the boundary value problems associated 
with a smooth elliptic operator in the region above a graph t = ip(x). When ip is C 1 , it 
was shown in [7] that the Dirichlet, regularity and Neumann boundary value problems 
are solvable with data in LP for 1 < p < oo, by the method of layer potentials. Our 
main theorem and the result of [5] show that the Dirichlet and regularity problems are 
solvable in this range of p when the boundary of the domain is defined by t = <p(x) 
where Vip G L°° f] VMO. The Carleson condition arises naturally here; the function 
!j = mo$ for the mapping $ : 1R™ — > {X — (x, t); t > 4>{x)} defined by (14.101) solves an 
elliptic equation in K™ with coefficients satisfying (II. ip . 

The paper is organized as follows. In Section 2, we give definitions and state our main 
results. Here we reduce the proofs of these main results to the case p = 2 and integer 
p for the regularity and Neumann problems, respectively. The proofs are applications 
of several general results, including the recent advances mentioned above. Section 3 
considers in detail the square function of the gradient of a solution. Section 4 explores 
the comparability of the square and non-tangential maximal functions and in Section 
5 the p = 2 version of our regularity result is established. In Section 6 we revisit 
bounds on the square function, this time by the co-normal derivative. Finally Section 
7 establishes a version of our Neumann result for integer values of p. 

Acknowledgements. J. Pipher was partially supported by NSF DMS grant 0901139. 
M. Dindos was partially supported by EPSRC EP/J017450/1 grant. D. Rule gratefully 
acknowledges the support of CANPDE. Part of this research was carried out during a 
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hospitality. 
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2. Definitions and Statements of Main Theorems 

Let us begin by introducing Carleson measures and the square function on domains 
which are locally given by the graph of a function. We shall assume that our domains 
are Lipschitz. 

Definition 2.1. Z C W l is an L-cylinder of diameter d if there exists a coordinate 
system (x, t) such that 

Z = {(x,t) : \x\ < d, -2Ld < t < 2Ld} 

and for s > 0, 

sZ:={(x,t) : |z| < sd, -2Ld < t < 2Ld}. 

Definition 2.2. Q C M. n is a Lipschitz domain with Lipschitz 'character' (L,N,Cq) if 
there exists a positive scale r and at most N L-cylinders {Z J }^ 1 of diameter d, with 
^ < d < Cor such that 

(i) 8Zj fl dfl is the graph of a Lipschitz function <pj, H^Hoo < L, <frj(0) = 0, 

(ii) on = U^ndn), 

j 

(Hi) Zj nn D Ux,t) e n : \x\<d, dist((x,t),dty < ^j. 

Definition 2.3. Let Q be a Lipschitz domain. For Q G dQ, X £ Q and r > we 

write: 

A r (Q) = dnnB r (Q), T(A r ) = nn b r (Q), 

5(X) = dist(X,dQ). 

Definition 2.4. Let T(A r ) be the Carleson region associated to a surface ball A r in 
dQ, as defined above. A measure \x inVt is Carleson if there exists a constant C = C(r ) 
such that for all r < r , 

MT(A r )) < Ca(A r ). 
The best possible C is the Carleson norm. When dfi is Carleson we write d[i G C. 
If lim C(ro) = 0, then we say that the measure fi satisfies the vanishing Carleson 

ro->0 

condition, and we denote this by writing d^i^Cy. 

Definition 2.5. A cone of aperture a is a non-tangential approach region for Q G dfl 
of the form 

F a (Q) = {X G Q : \X - Q\ < (1 + a) dist(X, dQ)}. 
Sometimes it will be necessary to truncate T a (Q), so we define T at h(Q) = T a (Q)nBh(Q) ■ 

Definition 2.6. If Q G R n , the square function of a function u defined on Q, relative 
to the family of cones {T a (Q)} Q( z d n, is 

S[a](u)(Q) = ( f \Vu(X)\ 2 (X)dist(X,d{l) 2 - n dx\ ' 
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at each Q G dfl. The non-tangential maximal function relative to {T a (Q)}Q e QQ is 

N a (u)(Q)= sup \u(X)\ 
xer a (Q) 

at each Q G dfl The truncation at height h of the non-tangential maximal function 
is defined by N[ a ^ h (u)(Q) = ^Pxer a (Q)nB h (Q) \ U (X)\, w ^h a similar notation Sja^/i for 
truncated square function. 

It will often be convenient to supress one or both of the parameters a and h in the 
square and non-tangential functions when their values do not play a significant role in 
an argument. So we may write S, S[ a ] or Sh to denote S[ a ] t h when no confusion should 
arise. Similarly we may abreviate Nt a ^h as N, N\ a ] or N^. 
We also define the following variant of the non-tangential maximal function: 

(2.1) N(u)(Q) = N [a] 4u)(Q) = sup (/ \u(Y)\ 2 dY 

xer a , h (Q) \Jb S(x)/2 (x) 

Definition 2.7. Let 1 < p < oo. The Dirichlet problem with data in L p (dfl,da) is 
solvable (abbreviated (D) p ) if for every f G C(dQ) the weak solution u to the problem 
Lu = with continuous boundary data f satisfies the estimate 

\\N(u)\\LP(dn,dcT) ^ \\f\\Lp(dn,da)- 

The implied constant depends only the operator L, p, and the Lipschitz character of the 
domain as measured by the triple (L,N,Cq) of Definition \2. H 

Definition 2.8. Let 1 < p < oo. The regularity problem with boundary data in 
if 1,p (<9f2) is solvable (abbreviated (R) p ), if for every f G H l > p (dVL) n C{dQ) the weak 
solution u to the problem 

Lu =0 in Q 
u\qb = f on dfl 

satisfies 

\\N(Vu)\\ LP(dn) < ||V t /||lp(«i). 

Again, the implied constant depends only the operator L, p, and the Lipschitz character 
of the domain. 

Definition 2.9. Let 1 < p < oo. The Neumann problem with boundary data in L p (dQ) 
is solvable (abbreviated (N) p ), if for every f G L p (dQ) D C(dQ) such that J dn fda = 
the weak solution u to the problem 

{Lu =0 in ft 

AVu -v = f on dfl 

satisfies 

\\N(Vu)\\ LP{dn) <\\f\\ LP{m) . 

Again, the implied constant depends only the operator L, p, and the Lipschitz character 
of the domain. Here v is the outer normal to the boundary dfl. 

We are now ready to formulate our main results. 
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Theorem 2.10. Let 1 < p < oo and let £1 c M be a Lipschitz domain with Lipschitz 
norm L on a smooth Riemannian manifold M and Lu = div(AVu) be an elliptic 
differential operator defined on Q with ellipticity constant A and coefficients which are 
such that 

(2.2) 6(X)~ 1 (osc B ( X ,s(x)/2)aij) 2 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(tq). Then there exists e = e(A,n,p) > such that if max{L, C(r )} < e then the 
(R) p regularity problem 

Lu =0 in Q 

u\an = f on dQ 
Jv(Vu) E L p (dtt) 

is solvable for all f with || Vr/||z,p(an) < oo. Moreover, there exists a constant C = 
C(A,n,a,p) > such that 

(2.3) \\N{Vu)\\ LPm <C\\V T f\\Lr(asi). 

In particular, if the domain Q is C 1 and A = (a^) satisfies the vanishing Carleson 
condition, then the regularity problem is solvable for all 1 < p < oo. More generally, if 
the boundary of the domain Q is given locally by a function <fi such that V0 belongs to 
L°° fl VMO, then, once again, the regularity problem is solvable for all 1 < p < oo. 

Proof. As follows from Theorem 15.21 the (R)2 regularity problem is solvable for opera- 
tors satisfying (12. 2p . provided e is sufficiently small. To complete the proof we will use 
[21 Theorem 1.1]. (There is also an older result by [T7] for symmetric operators which 
should be adaptable to the non-symmetric case.) 

According to [31 Theorem 1.1] (R)? solvability implies the solvability of (R)hs 1 (this 
is an end-point Hardy-Sobolev space boundary- value problem corresponding to p — 1). 
We also have from [31 Theorem 1.1] that under the assumption (R)hs 1 is solvable we 
have for p G (1, oo): 

(2.4) {R)p is solvable if and only if (D*) p > is solvable for p' = p/(p — 1). 

Here (D*) p / denotes the L p ' Dirichlet problem for the adjoint operator L*u = div(74*Vn). 

However by [51 Corollary 2.3] the L p ' Dirichlet problem for the operator L* is solvable 
under the assumptions of Theorem 12.101 (for sufficiently small e = e(p') > 0). Hence 
by (12 .4p the {R) p problem for the operator L is solvable proving our claim. □ 

As follows from the proof given above we also have a result for the endpoint p = 1. 

Corollary 2.11. Under the same assumptions as in Theorem \2.10\ the (R) H gi regular- 
ity problem for the operator L is solvable for all f with V^/ in the atomic Hardy space 
(c.f. [21 Theorem 2.3]j. Moreover, there exists a constant C = C(A,n,a) > such that 

(2.5) \\N(Vu)\\ LHd u)<C\\V T f\\n Hd uy 
For the Neumann problem we have the following: 
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Theorem 2.12. Let 1 < p < oo and let C M be a Lipschitz domain with Lipschitz 
norm L on a smooth Riemannian manifold M and Lu = div(AVu) be an elliptic 
differential operator defined on Q with ellipticity constant A and coefficients which are 
such that 

(2.6) ^{Xy 1 (osc B{x ,s(x)/2)aij) 2 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(r ). Then there exists e = e(A,n,p) > such that if max{L, C(r )} < e then the 
(N) p Neumann problem 

{Lu =0 inVt 

AVu -v = f ondQ 
N{Vu) E LP(dQ) 

is solvable for all f in L p (dQ) such that J dQ fda = 0. Moreover, there exists a constant 
C = C (A, n, a, p) > such that 

(2.7) \\N(Vu)\\ LP{dQ) <C\\f\\ LPm . 

In particular, if the domain Q is C 1 and A = (a^) satisfies the vanishing Carleson 
condition, then the Neumann problem is solvable for all 1 < p < oo. More generally, if 
the boundary of the domain Q is given locally by a function <fi such that V0 belongs to 
L°° n VMO, then, once again, the Neumann problem is solvable for all 1 < p < oo. 

Proof. As follows from Theorem 17.11 the ( N) p Neumann problem is solvable for opera- 
tors satisfying (17. ip . provided e is sufficiently small and p is an integer. To replace the 
condition (17. ip by (12. 6 p we use the same idea as Corollary 2.3] and Theorem 15.21 
For a matrix A satisfying (17. ip with ellipticity constant A one can find (by mollifying 
the coefficients of A) a new matrix A with same ellipticity constant A such that A 
satisfies (17. ip and 

(2.8) sup^X)- 1 !^ - A)(Y)\ 2 ; Y e B(X, 6{X)/2)} 

is a Carleson norm. Moreover, if the Carleson norm for matrix A is small (on balls of 
radius < r ), so are the Carleson norms of (17. ip for A and ( 12. 8p . Hence by Theorem 
17. II the {N)p regularity problem is solvable for the operator Lu = div(A'Vu). 

The solvability of the Neummann problem for perturbed operators satisfying (12. 8p 
has been studied in [13]. It follows by [131 Theorem 2.2] that the LP Neumann problem 
for the operator L is solvable, provided (12 .8p has small Carleson norm and the regularity 
(R) p and Neumann (N) p problems are solvable for L. Actually, the results in [13] are 
stated for symmetric operators, however careful study of the proof of [131 Theorem 2.2] 
reveals that symmetry is not necessary. 

However by Theorem 12.101 the {R) p regularity problem for L is solvable provided 
the Carleson norm of (17. ip is sufficiently small and {N) p Neumann problem for L is 
solvable by Theorem 17.11 Hence we have solvability of the Neumann problem (N) p for 
L by P3J Theorem 2.2]. 

If p > 1 is not an integer we use [121 Theorem 6.2]. (This result is also stated for 
symmetric operators, however, once again, symmetry is not necessary.) This theorem 
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implies that (N) p is solvable, provided (R)k and (N) k are solvable, where k is any 
integer larger than p. □ 

3. The Square Function for the Gradient of a Solution 

In this section we shall assume that Q is a smooth domain on a smooth compact 
Riemannian manifold M. As we shall see the case of a Lipschitz domain with small 
coefficients can be reduced to this situation via a pull back map of Dahlberg, Necas 
and Stein (see (OOjl ). 

We note that this case includes the most usual situation when Q C M. n is a bounded 
domain as in this case we can think of Q as being embedded into a large torus T n . 
We aim to establish local results near dQ. For this reason we introduce a convenient 
localization and parametrization of points near dfl. 

We want to write any point X G f2 near dQ as X = (x, t) where x G dQ and 
t > 0. The boundary dfl itself then will be the set {(y,0);y G dQ}. One way to get 
such a parametrization is to consider the inner normal N to the boundary dQ. The 
assumption that dQ is smooth implies smoothness of N. On Q we have a smooth 
underlying metric of the manifold M. 

We consider the geodesic flow Tt in this metric starting at any point x G dfl in the 
direction N(x). We assign to a point lefi coordinates (x, t) if X = TtX. This means 
that starting at x G dQ it takes time t for the flow to get to the point X. It's an easy 
exercise that the map (x,t) i— > X — Ttx is a smooth diffeomorphism for small t < to- 
Using this parametrization we consider the set Q to = {(x,t); (x,0) G dQ and < t < 
to}- 

Let us now deal with the issue of the metric. We want to work with the simplest 
possible metric on Q available. Since we only work on Q to we take our metric tensor 
there to be a product da ® dt where da is the original metric tensor on Q restricted to 
dQ. The product metric da (8) dt is different from the original metric on Q, but they 
are both smooth and comparable, that is the distances between points are comparable. 
Now we express the operator L in this new product metric. 

We note that under this pullback the new coefficients of our operator are going to 
satisfy the same Carleson condition as the original coefficients with Carleson norm 
comparable to the original. We observe in particular that the Carleson condition 
implies that V 'A G L^ c (Q to ) hence any solution of Lu = on Q to has a well defined 
pointwise gradient Vit. Furthermore, in the product metric da ® dt, the gradient Vw 
can be written as 

Vn = (V T u,d t u), 

where Vt is the gradient restricted to the n — 1 dimensional set dVt x {t = const}. 

Frequently throughout the paper it will be useful to localize to a single coordinate 
patch. The following definition gives a precise notion of coordinate frame. 

Definition 3.1. Let dQ be a smooth n—1 dimensional compact Riemannian manifold. 
We say that a finite collection of smooth vector fields (T T )™ =1 (m > n — 1) is an 
coordinate frame for dQ if: 
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• there is a finite collection of open sets U\, Ui, ■ ■ ■ , Uf. in W 1 1 and smooth diffeo- 
morphisms ipi : Ug — > dQ such that \J e fg(Ue) covers dVL, where Ug is an open 

subset of Ug such that Ug C Ug; 

• for each 1 < £ < k there exist a set Ag C {1, 2, ... , m} such that \Ag\ = n — 1 
and 

{^{f T )\~ reA e } = {£:\~,j = l,2,...,n- 1}. 

That is the pullback of the vectors T T to Ug, r G Ag restricted to Ug are just 
coordinate vector fields on Ug. 

Clearly, dQ has at least one such coordinate frame. Indeed, the existence of a finite 
collection (Ug, Ug, ipg) satisfying all assumptions of the previous definition follows from 
the fact that dQ is a smooth compact Riemannian manifold. Then on each Ug we 
consider vector fields ipgj-, j = 1, 2, . . . , n — 1 where ipg E C^(Ug) is a smooth cutoff 

function such that ipg\jj- = 1 and < ipg < 1 on Ug. Then 

{<pi*(ipg§:);l<£<k, l<j <n-l} 

is one such coordinate frame. Here tpg^ denotes the push-forward of a vector field from 
Ug onto dfl. 

We start with the following key lemma for the square function S(Vtu). 

Lemma 3.2. Let u be a solution to Lu = divAVu = 0, where L is an elliptic differ- 
ential operator defined on Q to with bounded coefficients such that 

(3.1) sup{5(X)|VMr)| 2 : Ye B 5(X )/2 (X)} 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(r ). 

Then there exists r\ > and K > depending only on the geometry of the domain 
Q, elliptic constant A and dimension n such that 

(3.2) f S 2 r/2 {V T u)do-~ [[ \V(V T u(X))\ 2 5(X)dX 

Jdn JJdQx(0,r/2) 

< K f \V T u\ 2 da + e [ N 2 (Vu)da + —\\Vu\\ 2 L2{n) , 
J aa Jan r 

for allr < min{ro, r\, to}- Here e > depends on the Carleson normC(ro) ande — > 0+ 

if C(ro) — > 0. iV r denotes is the non-tangential maximal function truncated at height 

r, S(X) = t for a point X = (x,t) G dVt x (0,to)? Vt« is the tangential gradient of u 

on dfl x {t} and dX is the product measure dadt. 

Proof. In order to establish ( 13. 21) we localize to coordinates. Let U = Ug be one of the 
sets from Definition 13.11 with corresponding map ip = ipg, equally set U = Ug. We can 
now consider the operator L as being defined on an open subset U x (0, to) of , where 
dQ corresponds to the hyperplane {(x, 0); x G U}. We achieve this by pulling back the 
coefficients of L from Q ta to U x (0, t ) using the smooth map $ : (x, t) \-> (<p(x),t). At 
this stage we also pull back the product metric da (g> dt from dQ x (0, to) to U x (0, to) 
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and we get another product metric that we (in a slight abuse of notation) still denote 
by da ® dt on U x (0, t ). 

Since we are going to use a partition of unity we also consider a smooth cutoff 
function (f)(x,t) = (j){x) defined on M™, independent of the t = x n variable such that 

< <f>(x, t) < 1, supp <j> C U x R. 

Instead of the left-hand side of ( 13. 2ft (by the ellipticity of the coefficients) we are 
going to estimate a similar object 

(3.3) ff ^-{diVXdjV^tdadt, 

JjUx(0,r) a nn 

for functions V = Vu ■ T T , for 1 < r < m. Here and below we use the summation 
convention and consider the variable t to be the n-th variable. We begin by integrating 
by parts 



(3.4) ff ^(d l V)(d j V)(j)tdadt = -f d (\V\ 2 )^<ptUida- 

JjUx(0,r) a nn 2 Jt/ X { r } a nn 

- ff —V{LV)ct)tdadt-ff V{d j V)a ij d i (—\ dadt. 

JJux(0,r) a nn JJux(0,r) \ a nn J 

Here v$ is the z-th component of the outer normal is, which (given we consider a product 
metric) is is just the vector e n for the boundary U x {r}. Hence the first term is non- 
vanishing only for i = n. We work on the last term, as it is the most complicated. This 
one splits into three new terms, one when the derivative hits t (where only the term 
with i = n will remain) and another two when it hits (ft and l/a nn : 

~ ff V(djV)^<j)dadt - ff V(d,V)^(d l( j))tdadt 

JJux(0,r) a nn JJux(0,r) a nn 

(3.5) + ff Vid^^fidtdnn^tdadt. 

JjUx(0,r) a nn 

Consider now the first term of (13.51) . For j — n, as is independent of x n = t, we only 
get 

(3-6) -Iff d n (\V\ 2 <P)dadt= l - f \V\ 2 <pda- 1 - f |Vf0da 

2 JJux(0,r) Z JU Z JUx{r} 
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For j < n the first term of (13. 5p is handled as follows. We introduce an artificial 
term 1 = d n t inside the integral and integrate by parts. 



" Iff d ] {\V\ 2 )^ L <P(d n t)dadt = -\( d 3 {\V\ 2 )^da 

z JJUx(0,r) a nn & Jljx{r} a nn 

(3-7) + \ff d n (d 3 (\V\ 2 )—A tdadt = -lf d 3 {\V\ 2 )—4>tda 

* JJUx(0,r) V a nn / * JUx{r} a nn 

+ l [[ dM\V\ 2 )—^dadt+ l -(f d 3 (\V\ 2 )d n (^A^tdadt. 

* JJUx(0,r) a nn * JjUx(0,r) \ a nn J 

The first term here gets completely cancelled out by the first term of f l3.4[) as they have 
opposite signs. The second term can be further integrated by parts and we obtain 

Iff dM\y\ 2 )—4>tdadt = -Iff WVfidjf^A&dadt 

^ JJux(0,r) a nn * JJux(0,r) \ a nn / 



1 ff i v 



(3.8) - - // d n {\V\ 2 ) — {dj<t>)tdadt 

JJUx(0,r) a nn 

We now notice that the last term of (I3.5p . the third term on the righthand side of 
(13. 7p and the first on the righthand side of (13. 8p are of same type and are bounded 
from above by 



(3.9) C // \V\\VV\\VA\(f)tdadt. 

Here VA stands generically for either Va nj -, Va nn . Estimating ( 13. 9p further we see 
that, using the Cauchy-Schwarz inequality, ( 13.91) is less than 

(3.10) c(ff \V\ 2 \VA\ 2 (j)tdadt] ( ff \VV\ 2 (j)tdadtY . 

Using the Carleson condition on the coefficients, and the fact that the Carleson norm 
is less than e we get that this can be further written as 

(3.11) Ce 1/2 ( f N r (V) 2 dy] ' ( ff \VV\ 2 (j)t dadt] ' 
< — ff \W\ 2 <j)tdadt + — f N r (V) 2 dy 



z JJUx{0,r) * JU 

where the last line follows from the inequality between arithmetic and geometric means. 
We observe that the first term on the second line is no more than one half of (13.31) and 
hence can be incorporated there. 

Let us summarize what we have. For some constant C > and e (the Carleson 
norm) we have that 
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(3.12) If ^-(diVXdjV^tdadt 

JjUx(0,r) a nn 

< Ce [ N r (Vfdy+ [ \V\ 2 0da- [ \V\ 2 <pda + 

JU JU JUx{r} 

+ [ d n {\V\ 2 )<f>tda- [[ —V(LV)0tdadt + E. 

JUx{r} JJux(0,r) a nn 

The fourth term on the righthand side is the first term of (13.41) for i = j = n. Here 

(3.13) E = -[[ d 3 (\V\ 2 )^(d4)tdadt- [[ d n (\V\ 2 )^(d j( f))tdadt. 

JJux(0,r) a nn JJux(0,r) a nn 

We call E "the error terms" these are the second term of (13.51) and the second term 
on the righthand side of (13. 8p . Both terms are of same type and contain <9;0 for i < n. 
(Recall that d n <p = 0). 

At this point we have to use the fact that V = V« • T i; for 1 < i < m, where (T i )™ =l 
is a frame from Definition 13.11 It follows that in our local coordinates 

V = '^^b k Vk, for some smooth functions b k on U. 

k<n 

Here 

du 

Vk = d k u = - — , for k = 1, 2, . . . , n - 1. 

ox k 

We denote by v n = d t u. We observe that each v k is a solution of the following auxiliary 
inhomogeneous equation: 

(3.14) div(AVv k ) = Lv k = -div((«9,A)v) = divF fc , 

where the i-th component of the vector F k is (F k ) 1 = —(d k a,ij)djU = —(d k ciij)vj. 
It remains to deal with the second term of the last line in (I3.12p . Clearly, 

(3.15) LV = [<9iK(^& fe )H + a^djb^diVk + a^ill/ U),,;, + b k Lv k ] . 

k<n 

We will have to deal with these four terms. We start with the second and third 
ones as they are the easiest. We observe that since b l are smooth, both b % and V6 l 
actually satisfy the vanishing Carleson condition. Hence these two terms put into the 
expression 

(3.16) [[ —V(LV)<ptdadt 

JjUx(0,r) a nn 

can be estimated by 

(3.17) CV [[ \Vu\\Vv k \\B\4>tdadt, 

k<n JJux(o,r) 
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where B stands for a generic coefficient such as aij(djb k ) or ciij(dib k ) with \B\H being 
a density of a vanishing Carleson measure. Hence in the same spirit as we dealt with 
(EU) we get 



(3.18) V // \Vu\\Vv k \\B\<ptdadt< 

k<n JJux(0,r) 

< KJ2 ff \Vv k \ 2 <ptdadt + C{K)e f N r (Vu) 2 dy. 

k<n JJux(0,r) JU 

The K is this formula can be arbitrary small, e only depends on the Carleson norm of 
\B\H near the boundary, hence it can also be arbitrary small by making r\ > smaller 
if necessary. We choose K sufficiently small so that the first term on the second line 
of (I3.18P can be hidden on the left-hand side of (13. 2p . 

Next we look at the first term of (I3.15P as we place it into (I3.16p . We obtain 



(3.19) 



ff —Vd i (a ij (d j b/'))v k <frtd(Tdt<cff |Vm| 2 |V 'B\<f>tda dt. 

JjUx(0,r) a nn JJux(O.r) 



'Ux(0,r) 



Here |V-B| 2 t satisfies the Carleson condition with small norm since |Vajj| 2 £ has a small 
Carleson norm, and djb k and its higher derivatives are smooth functions (and so they 
satisfy the vanishing Carleson condition). Hence by the Cauchy-Schwarz 

1/2 / rr \ 1/2 



\Vu\ 2 \VB\<f)tdadt < ( ff \Vu\Hdadt] iff \Vu\ 2 \VB\ 2 tdadt] 

\JJux(0,r) J \JJux(0,r) J 



Ux(0,r) \JJUx(0,r) J \JJUx(0,r) 

(3.20) < Cre l/2 f N 2 (Vu)da. 

Ju 

Here we observe that the last term on the first line is of the same type as the first term 
in ( 13.1 Op we have handled before. It follows that this term is small even if the Carleson 
condition has a large norm by choosing r small. 

We now handle the last term of (13. 15[) using ( I3.14p . Placing this into (I3.16P yields a 
term 

(3.21) ff —Vb k d i {{d k a lj )v j )(j)td(jdt = 

k<n JJux(0,r) a nn 

rr fb k (j) \ f b k d k CL 

= — 2_. // 9% ( Vt) (d k ciij)vj dadt + / —Vvj<f)tdcr, 

k<n JjUx(p,r) \ a nn J JUx{r} a nn 
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where we integrate by parts and only obtain a boundary term when % = n. We deal 
now with the solid integral. This gives 

(3.22) -Jiff dJ—^idka^Vv^tdadt- 

k<n J->Ux(0,r) \ a nnJ 

-JZff b -^^d t Vv^tdadt- 

k<n JJux(0,r) a nn 

- V ff b !^lVv J (d t( f ) )tdadt - 

k<n "ux(0,r) a nn 
k<n JJux{0,r) a nn 

where the last term only appears for % = n (as d n (t) = 1). We notice that the 
first term here is of the same type as the first term in ( 13.101) and hence bounded 
by e fjj N r (Vu) 2 da. The second term is handled exactly as (I3.9P (noticing that \vj\ < 
|Vw|). Hence this term is (in absolute value) no greater than 

K ff \VV\ 2 (j)tdadt + C{K)e f N?(Vu)dy, 

where K > can be arbitrary small and e > depends on the Carleson norm of the 
coefficients. Thus as before by choosing K sufficiently small this term can be absorbed 
into the left-hand side of ( I3.12p . 

The third term of (I3.22p is another "error" term of type similar to ( I3.13p . We will 
handle this at the end. Hence the only term remaining is 

_ ^ ff ^m Vvj f dadt = - W/ b -^Vv j cf>d n (t) dadt. 

Here we introduced an extra term 1 = d n (t) and now integrate by parts again. This 
gives 

(3.23) ~ Hf h ^^Vv 3 <\>tdodt + 

k<n Jux{r} «nn 

+ ff d n ( ——) (d k a n j)Vvj 4>t dadt + 
+ J2 ff b -^^ivd n v j <j>t dadt + 

k<n "ux(0,r) a nn 



+ V ff h ^^i dnVv . t dadt + 

'Ux(0,r) a nn 



k<n 



+ V // bk{9Aa ^ Vv j( j>tdadt. 

k<n JJux(0,r) a nn 



14 



MARTIN DINDOS, JILL PIPHER AND DAVID RULE 



The first four terms are of same type we have encountered before. The first term here 
is cancelled by the last term of (I3.2ip . The second term is bounded by e J v N 2 (Vu) da 
(c.f. (I3.10p ). The third term is like flUTTJ) and the fourth like the second term of (13.221) . 
Finally, in the last term we have two derivatives on the coefficient (d n dka n j) but only 
one of the derivatives is in the normal direction since k < n. Hence we integrate by 
parts one more time (moving the d k derivative). We get 

nd n d k a nj ) Vv ^ tdadt= 

'Cx(0,r) a nn 



(3.24) 



E 

k<n 



E 



k<n 



Ux(0,r) V^nn 

b k d v a. 



(d n a n j)Vvj (fit dadt 



J n KJj nj 



k<n JJUx (°< r } Unn 

E 



E 



k<n JJ Ux(0,r) u nn 

b k d„a. 



VdkVj (fit dadt 



dkVvj (fit dadt 



Vvj (d k (fi)t dadt. 



k<n uJUx (°> r ') nn 

Here the second, third and fourth terms are like the second, third and fourth terms in 
(I3.23P and are handled likewise. Finally, the last term is another of the "error terms" . 
This concludes the analysis of the term (I3.16P in (I3.12p . 

Finally, we sum over all choices of functions V = V T = Vw • T T , for 1 < r < m and 
over all sets Ui (from Definition 13. 1 1) choosing the smooth cutoff functions (fi = (fig in 
( 13. 3p such that they are the partition of unity, that is 



E' 



1 on dQ 



and 



supp 



cu t . 



We first observe that the terms we called "error terms" completely cancel out. There 
are the terms in (13.131) plus two extra terms later on. This is due to the fact that 



0. That means that summing over r these terms equal to zero. This 



cancellation happens even if we work on different coordinate charts since the term we 
started our calculation (13.31) does not depend on choice of coordinates. Hence after 
taking into account all remaining terms we have by (13 . 1 2p : 



(3.25) 



\V(V T u(X))\ 2 S(X)dX 



dflx(0,r) 



« E // \VV T (X)\ 2 5(X)dX< 

i=1 JJdnx(0,r) 

< I<[ \V T u\ 2 da + e [ N 2 {Vu)da 
Jdn Jan 



'an 

m r 



+ A 'E 



T=l 



on 

2\ 



d n {\V T \y da - / \V T \ 2 da 

dnx{r} Jdnx{r} 
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At this point we have to deal with the last two terms 



f d n {\V T \ 2 )r da - f \V T \ 2 da = 

JdQx{r} JdQx{r} 

= f d n {\V T \H) da - 2 ( \V T \ 2 da< f d n {\V T \ 2 t) da. 

Jdnx{r} JdQx{r} J dflx{r} 

We would like to estimate this by a solid integral by integrating r over an interval 
(0, r') and averaging. This yields 

r < 

1 ' ' " _ / IT/ |2 



(3.26) - / / d n {\V T \H)dX= / \V T \ 2 da. 

r J0 Jdnx{r} Jdflx{r'} 

This term is still not a solid integral so we use the averaging technique one more 
time by integrating over r' and averaging over an interval (0,ro). This yields a solid 
integral 

- ff \V T \ 2 dX<-[[ \Vu\ 2 dX. 

Going back to (13.251) we have to perform this double averaging procedure on all terms. 
This leads to introduction of some harmless weight terms and finally an estimate 



(3.27) // \V(V T u(X))\ 2 6(X)dX < 

< K I \V T u\ 2 da + e( N 2 {Vu)da + — 1| Vw||| 2(n) . 
Jan Jan r 

□ 

Lemma 13.21 deals with square function estimates for tangential directions. We have 
following for the normal derivative: 

Lemma 3.3. Let u be a solution to Lu = divAVu = 0, where L is an elliptic dif- 
ferential operator defined on Q to with bounded coefficients which are such that is 
the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(r ). Then 



(3.28) f S 2 (d n u)da= ff \V(d n u(X))\ 2 S(X)dX < 

Jan JJanx(o,r) 

< K ff \V{V T u{X))\ 2 5{X)dX + e f N 2 {Vu)da 
JJanx(o,r) Jan 

= K f S 2 .{V T u)da + e f N 2 {Vu)da 
Jan Jan 

provided r < min{ro,to}- Here e > depends only on the Carleson norm C(ro) and 

e — > 0+ as C(ro) — > and K only depends on the domain, ellipticity constant and 

dimension n. 
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Proof. We integrate by parts in dQ x (0,r). We use the notation introduced above 
where we denoted v n = d n u. Clearly 

(3.29) ff \Vv n {X)\ 2 5{X)dX 

= ff \V T v n (X)\ 2 5(X)dX+ [[ \d n v n (X)\ 2 5(X)dX = 
JJanx{o,r) JJanx(o,r) 

= ff \d n (y T u(x))\ 2 5{x)dx+ ff \d n v n {x)\H{x)dx 

JJanx{o,r) JJanx(o,r) 

The first term is clearly controlled by the square function of Vtu. It remains to deal 
with the second term. Since 

\a n nd n v n \ 2 = \d n (a nn v n ) - d n (a nn )v n \ 2 < 2\d n (a nn v n )\ 2 + 2\d n (a nn )v n \ 2 . 

We see that by the ellipticity assumption 

(3.30) ff \d n v n (X)\ 2 5(X)dXn ff (a nn (X)) 2 \d n v n (X)\ 2 5(X)dX < 

JJdnx(0,r) JJdnx(0,r) 



< 2 \d n (a nn v n )\ 2 tdX + 2 \d n (a nn )v n \ 2 tdX. 

JJdQx(0,r) JJdnx(0,r) 

Here as before X = (x,t), i.e. t is the last n-th coordinate. The second term (using 
the Carleson condition) is bounded by e J dn N 2 (Vu) da. We further estimate the first 
term. Using the equation u satisfies we see that 

d n (a nn Vn) = - ^2 diiOijOjU). 
From this point on we use local coordinates. It follows that 



(3.31) ff 
JJa 



\d n (a nn v n )\ 2 tdX < (n 2 -l) V // ^{a^u^H dX 



2(n 2 -l) V ff \d i {a lj )\ 2 \d J u\HdX+ ff \a id | 2 |^-M| 2 t dX 

{i,j)^(n,n) \-"dVx(0,r) JJdnx(0,r) 



The first term here is of the same type as the last term of f 1 3 . 3 j) and is bounded by 
e J dn N 2 (Vu) da. Because (i,j) ^ (n,n) 

l&djul 2 < |V(V t m)| 2 , 

hence the last term of (I3.3ip is also bounded by the square function of Vt«- D 

4. Comparability of the Nontangential Maximal Function and the 

Square Function 

If we combine the results of Lemma 13.21 and 13.31 we obtain the following inequality. 

Lemma 4.1. Let u be a solution to Lu = div AVu = 0, where L is an elliptic differen- 
tial operator defined on Q to with bounded coefficients which are such that Ii3. 1\) is the 
density of a Carleson measure on all Carleson boxes of size at most ro . 
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Then there exists T\ > and K > depending only on the geometry of the domain 
Q, ellipticity constant A, dimension n and the Carleson norm of coefficients such that 

(4.1) / S 2 r/2 {Vu)da <K [ N 2 (Vu)da, 

Jon Jan 

for all r < min{ro, rr, to}- 

Proof. We observe that first two terms on the righthand side of (13.21) can both be 
bounded by K f Qn N 2 (Vu)da. Recall that the last term — ||Vw||| 2 (n) appears there 
due to averaging of (I3.26p . This last averaging is however unnecessary as the righthand 
side of (13.261) can be directly bounded by a multiple of J dfl N 2 (Vu)da. From this (14. ip 
follows. □ 

We note that Lemma [4.11 holds regardless of the size of Carleson norm C(r ). We 
would like to establish an analogue of Lemma |4~T1 for values p different from 2. In order 
to do that we first observe that a local version of Lemma [4. II is also true: 

Lemma 4.2. Consider an operator L defined on a subset 2U x (0,r) of R" , with 
r ~ diam(U). Then there exists K > depending only on the ellipticity constant A, 
dimension n and the Carleson norm of coefficients such that 

(4.2) f \V 2 u\tdadt<K [ N 2 {Vu)da. 

JUx(0,r) J2U 

Proof. The proof is essentially same as the proof of Lemma 14.11 since the estimate 
(13.21) is based on local considerations. However, the terms of type (13.131) have to be 
considered now as they only disappear in the global estimate. Observe that \di<p\ < C/r 
hence these "error" terms are bounded from above by 

C [[ \V 2 u\\Vu\±dadt. 

JJux(0,r) 

By Cauchy-Schwartz this can be further bounded by 

/ rr \ V 2 / rr \ 1 / 2 

C[ // \V 2 u\ 2 tdadt) // \Vu\ 2 ±dadt 

Since \Vu(X)\ < N(Vu)(Q) for all X e T(Q) the term JJ Ux(0r) \S7u\ 2 ±da dt is further 
bounded by 

From this (S2D follows. □ 

We claim that Lemma 14.21 implies that the square function is controlled by the non- 
tangential maximal function in LP for p > 2 as well. 

Lemma 4.3. Let u be a solution to Lu = div AVu = 0, where L is an elliptic differen- 
tial operator defined on Q to with bounded coefficients which are such that Ii3.1\) is the 
density of a Carleson measure on all Carleson boxes of size at most r . 
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Then for any p > 2 there exists r\ > and K > depending only on the geometry of 
the domain Q, ellipticity constant A, dimension n and the Carleson norm of coefficients 
such that 

(4.3) f S p r JVu)da<K f N^{Vu)da, 
for all r < min{r , r 1; t Q }. 

Proof. The lemma has already been proved when p = 2, since then it is just the 
statement of Lemma 14. 1\ so we only need to consider p > 2. Moreover, it suffices to 
prove (I4.3P on each coordinate patch Ut for I — 1, 2, . . . , k. In fact, we can go slightly 
further and say it is sufficient to prove 

(4.4) f S p r JVu)da<K [ iV r p (Vu)Ax 

Ju° Ju° 

for each £, where Ut C C Ut- Because we only need to consider p > 2, Lemma 
2 on page 152 of [18] shows that to prove (I4.4p it is sufficient to show the relative 
distributional inequality 

\{xe U° | S Hr/2 (Vu)(x) > 2A, M(N r (Vu) 2 )(x)^ < a \}\ 

< Ca 2 \{x G U? | S [2a]>r/2 (Vu)(x) > X}\, 

where M is the Hardy-Littlewood maximal function on U®. 

Now, for each £, we describe a localised Whitney decomposition of the set where 
5"[2a],r/2(Vn) > A (c.f. [HI A-34] , which we follow here). First for each £ we find a 

finite number of cubes Qtj for which Ut C UjQt,j Q Ut and the side length £(Qt,j) 
of Q^j is comparable with r. We denote Ptj = {x G Q^j | 5 , [ 2 a],r/2(VM)(x) > A} and 
Ktj = {x G | ^.^(V^i) < A} 

Fix a pair (£, j). If K £ j is empty, define Ttj := {Q^j}- If -K^j is non-empty we will 
define Ttj to be a collection of dyadic sub-cubes of Qtj in the following way. First 
observe that we can write P £ j as the union of 

P£. = {x G P Aj | 2£(g^)v^2- fc < dist(x, A^) < M(Qt tJ )VE2~ k } 

for k eN. 

We can find 2 n_1 dyadic sub-cubes of Qtj by bisecting each side of Qt,j- We denote 
the collection of these 2 n_1 cubes as D\- and each cube in the collection has side length 
equal to £(Qej)/2. Equally, we can find 2 n_1 dyadic sub-cubes of each cube in D\- 
by again bisecting each side of it. Thus, we have 2 2 ( n_1 ) subcubes of the cubes in 
Djj which have £(Qtj)/2 2 . We denote the collection of these 2 2(jl ~ l ' ) cubes by D\ -. 
Continuing inductively D k j is a collection of 2 fc ( n_1 ) dyadic cubes with side length equal 
to £(Qt,j)/2 k . 

Let T[ j be the collection of all cubes Q in D k - for some keN such that QnPfa ^ 0. 
Let Q G T't j and pick x G Q fl Pjfj. Observe that 

£(Qtj)2- k y/n~^l = dist(x, K ej ) - £(Qij)2- k Vn~^l = dist(x, K e ,j) - £(Q)V^1 
< dist{Q,K t>j ) < dist(x,K £J ) < 4£{Q^)2- k Vn - 1 
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and so, 

(4.6) £(Q)Vn^T < dist(Q, K e ,j) < M(Q)V^T. 

Given that J 7 ^- is a collection of dyadic cubes, any two cubes which intersect have the 
property that one is contained in the other. Thus, we may define Ttj to be the set of 
cubes Q G T[j such that if Q' G T[ j and Q PI Q' ^ 0, then Q' C Q. That is Tt# is the 
set of maximal cubes in T[ y Clearly then, Ttj is a collection of disjoint dyadic cubes. 
Our Whitney decomposition of {x G Ug | 5[ a ] )T ./ 2 (Vu)(x) > A} is then the collection 

Ti> := \jTt,y 

This collection has the properties that each Q G J-} is such that either (14. 6 p holds or 

£(Q) ~ r, 

|J Q = {x G UjQ^j I S[2a],r/2(Vu)(x) > A}, 
QeT t 

and there exists a constant C such that there are at most C cubes that intersect at 
any given point. 
Fix Q G Ti and set 

R := {x G Q | S r/2 (Vu)(x) > 2A,M(iV r (Vw) 2 )(x)2 < «A} 

If x G -R and ( 14. 6 p holds for Q, then there exists x' such that dist(x,x') < 4y/n£(Q) 

and Spa] r/2(Vu)(x') < A. Consequently there exists a constant a such that 

(4.7) 

S [aUi{Q) (Vu) 2 (x) > S Hr/2 (Vu) 2 (x) - [[ \V 2 u\H 2 - n dadt 

i/r [a]>r , /2 (a;)n(R»-ix(a^(Q),r/2)) 

> S [a]: r/ 2 (Vu) 2 (x) - S [2a ],r/2(^u) 2 (x') 

> 4A 2 - A 2 = 3A 2 

Then, if R is non-empty (say xq G R), we can apply Lemma [4.21 to conclude that 
\R\ < ^ / S [aW( Q)(Vw) 2 da < — / |V 2 w| 2 tcM 

<^T I N am {Vufda < 2nC * lQl M(N r (Vu) 2 )(x ) < TCKa 2 \Q\. 

Furthermore, if ( 14. 6p does not hold, then r/2 ~ ^(Q), so we may repeat ( 14. 8ft without 
the need for (14.71) . Finally, we observe that the inequality \R\ < 2 n CKa 2 \Q\ is trivial 
if R is empty. Thus, summing over Q G Tn we obtain (I4.5p with [7° = UjQej. □ 

Now we would like to establish a reverse inequality showing that the non-tangential 
maximal function can be dominated by the square function. As we shall see in the 
proof we will have to assume small Carleson norm. We start with the following local 
lemma working in coordinates on with boundary IR™" 1 . 

Lemma 4.4. Let u be a solution to Lu = divAVu = 0, where L is an elliptic differen- 
tial operator defined on Wl with bounded coefficients such that A3. 1\) is the density of 
a Carleson measure on all Carleson boxes of size at most r with norm e. 
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Let (J) be a non-negative Lipschitz function and let Q be a cube in R™" 1 with r = 
diam(Q). Suppose that 0(x) < 12r /a for x G Q* . Here Q* is a dilated Q by factor 
of 5. Then if \\V4>\\l^(r»-i) is sufficiently small, there are constants a (c.f. Definition 
\2. 5\) and C (depending only on A, Carleson norm e, || V</>||£oo( R n-i) and a), such that 

l|Vu(.,#.))|& w) < C (||5(V«)||i a(g . ) +e||JV(V«)||i aw . ) 

(4.9) + ||N(V«)|| ii! i W . ) ||5'(V«)||x.(g.) + r B - 1 |V«(X P )| a ), 

where X r is an arbitrary corkscrew point, i.e., any point in {X = (x,t); <ft(x)+r/2 <t< 
4>(x) +Qr/a}. The square and non-tangential maximal function in \4-9$ are defined 
using non-tangential cones r a (.). Both square function and non-tangential maximal 
functions on the righthand side can be truncated at a height that is a multiple of r. 

Proof. Recall the mapping $ : R™ — > VL^ = {X = (x,t);t > 4>{x)} used by Dahlberg, 
Keing and Stein (see for example [2] or [15] and many others) defined as 

(4.10) <5>(X) = (x,c t+(6 t *<p)(x)), 

where {6 t )t>o is smooth compactly supported approximate identity and Co can be chosen 
large enough (depending only on || V0||£°°( R n-i) so that <3> is one to one. We pull back 
the solution u in Vt^ of div(A Vu) = to a solution v = u o $ of a different second order 
elliptic equation div(B'Vv) = 0. 

The coefficient matrix B satisfies ellipticity condition with constant that is a multiple 
of A and which depends on ||V0||ioom»-i). Also if e is the Carleson norm of 

sup{t|Va y (F)| 2 : Y 6 B t/2 ((x, *))}, 
then the Carleson norm of 

sup{t|V%(r)| 2 = Y e Bt /2 ((x,t))}, 
for B = (bij) will only depend on e and ||V0||loo. Furthermore, if ||V0||£°o is small 
enough, then the Carleson norm of the matrix B can be guaranteed to be less than 2e. 

We choose a smooth function £i : W 1 ^ 1 — > R such that £,i(x) = 1 for x G Q, 
Wi\ ^ 16/r and support contained in a concentric dilation (9/8)Q. Choose another 
function £ 2 : [0, oo) — > R such that ^ 2 (t) = 1 on [0, r], \£' 2 \ < 5/r and support contained 
in [0,2r]. Now define £(X) = £{x,t) = 

Denote by Wi = diV for % — 1, 2, . . . , n. For each % < n — 1 we have 

f Wi (x^fUx)dx = - [[ d n {wH){X)dX 

(4.11) = - // 2w t (d n w^dX- [[ wfa&dX. 

The second term on the right-hand side of (14. lip is controlled by r" 1 JJ K wf where 
K = {X = (x,t);x G Q* and r/3 < t < 7r/a}. Let X r be any point in K and choose 
K' and K" to be the appropriate concentric enlargements of K. We set c = ff K , Wi. 
Using [SI Thm 8.17] we may further estimate this term by 
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r 1 jj ( Wi - Wi{X r )f dX + r" 1 jj w 2 (X r )dX 

< Cr 11 - 1 osc K (wi) 2 + CV^MX,.)! 2 

< Cr™" 1 sup \ Wi - c\ 2 + Cr"" 1 \wi(X r ) | 2 



A" 



< Cr' 1 ^ | Wj - c| 2 dX + Cr n - 1+2 ( 1 - n / 9 )||(9 i 5)w|| 2 q(A -, ) + Cr 71 " V^)] 2 , 

for q > n. Here we are using (I3.14p with matrix A replaced by B. Using the Poincare's 
inequality and Carleson condition for B this can be further estimated by 



C // iViw^rdX + Ce^NiVu^l^ + Cr^WiiXr)] 2 

JJ K" 

(4.12) < C|| < S(V U )||| w + C £ 2 ||7V(V W )||i 2(Q , ) + Cr"-> i (X r )| 2 . 
The first term on the righthand side of (14. lip can be estimated by 

- ff 2w i {d n w i )£dX 
JJr™ 

= - ff 2w i (d n w i )£d n (t)dX = 2 ff [d n (wi(d n Wi)£)]tdX 

jjr" jjr^ 



= 2// {dnWiYZtdX + 2 w i (d 2 n w i )&dX + 2 Wi {d nWi )^ 2 tdX 
JJr™ JJr™ JJr™ 

=: I + 11 + III. 

Using the fact that i < n — 1 we see that d 2 Wi in the term II can be written as did n w n . 
This gives 



II = -2 // (d n w n )di(wiZ)tdX 
JJr™ 

-2 ff {d n w n )d i {w i )itdX-2 ff (d n w n )wA(Zi)&dX 
JJr™ JJr™ 



Hi + H 2 . 
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We observe that the terms I and Hi are both bounded by the square function 
||'S'2r( w )|li2(Q*)- This is further bounded by ||S , (Vm)||^ 2 (q*), where the square func- 
tion is truncated at a greater height or not truncated at all. For II2 and III we have 
C 



II, + III < 



I Vw w t dX 



*x(0,2r) 



< C t 

"x(0,2r) 

< C\\S(w)\\ L 2 m 



WwrtdX 



1/2 



w 



1 

, r 



*x(0,2r) 
2r \ 1/2 

Iwl 2 dt dx 



1/2 



< C\\S(w)\\ L * {Q , 



— \N(w)\ 2 dx 



1/2 



C\\S(w)\\ L2m \\N(w)\\ L * m . 



This bounds (14. lip by terms that appear on the righthand side of (14. 9p . 

It remains to estimate L„_i w n (x, 0) 2 £,i(x)dx. We estimate instead an expression for 
co-normal derivative H = 52 . b n jWj. This is sufficient since 



(4.13) 



< n 



%(i,0) 2 ^i(i)cfiRi / (b nn w n ) 2 (x,0)£ 1 (x)dx 
/ H 2 ^i dx + / (b n jWj) 2 !;i dx 

< n f H 2 ^dx + CY] f w 2 (x,0)Ci{x)dx 

Hence if we can obtain estimates for the first term we are done since the second term 
has already been bounded. We proceed as before. 



H(x,0) 2 ^(x)dx 



(4.14) 



ff d n (H 2 0(X)dX 
J Jr.™ 

ff 2H(d n H)£dX- ff H 2 ^' 2 dX. 
JJr™ JJr™ 



As before we observe that the second term can be bounded by r~ l JJ K w 2 . The 
calculation we have done above holds for any i even i = n giving us bound (14.121) . 
It remains to deal with the first term. Using the equation div(-BV-u) = 

d n H = ^ d n {b nj djv) = ~y~] diibijdjv) = d^b^Wj). 
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It follows that 



- ff 2H(d n H)£dX 
JJwq_ 

(4.15) = ff 2Hd i (b ij w j )(d n t)^dX = -J2 ff 2d n {Hd i {b ij w j )OtdX 

= -V ff 2{d n H)d i {b ij w j )it dX - ff 2Hd i d n {b ij w j )£tdX 
%<n JJv» JJr- 

- ff 2Hd i (b ij w j )£ 1 &t dX = l + Tl + ILL 

As before we do further integration by parts for the term II. 
S = ff 2d n {b i fw j )d i {H£)tdX 



= 2 ff 2d n (b ijWj )(d i H)&dX+ ff 2d n (b ij w j )H(d i ^tdX 
JJr™ JJr™ 

= Ui+fh. 

We observe that when the derivative in terms II 2 and III does not hit the coefficients 
bij these can be estimated exactly as the corresponding terms II2 and III. When the 
derivative falls on the coefficient we get "error terms" that can be estimated using the 
Carleson measure property of the coefficients. In particular the term from III is of the 
same form as (13.201) and is handled analogously. The term we obtain from II 2 is of a 
different nature and can be bounded above by 

\w\ 2 \VB\±dX. 

r Q*x(0,2r) 

By Cauchy-Schwarz this is no more than 

C ( ff \VB\ 2 \w\Hdx) V2 ( ff |w| 2 ^ dx) V2 

< Ce\\N(w)\\ L , m (J l -J* r \w\ 2 dtdi\ 

( f 2r \ 1/2 

< Ce\\N(w)\\ L * m ^ -\N(w)\ 2 dx) = CeH^w)!^^. 

The terms I and Hi contain both a derivative acting on H and a derivative acting on 
bijWj. We deal with these in two parts: (a) when the derivative acting on H = b n jWj 
falls on b n j and (b) when it falls on Wj. First we deal with case (b). When the derivative 
acting on b^Wj does not hit the coefficients, we can handle them as the corresponding 
terms I and Hi. When this derivative falls on the coefficients, the term we get from I 
is again of the same nature as (I3.20|) and the term we get from Hi looks like (13.91) . so 
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these terms are handled as before. Finally we deal with ), where we either get 

terms of the form (13. 9p . which we have dealt with before, or terms of the form 

(4.16) ff \w\ 2 \VB\ 2 t£dX <e [ N 2 {Vu)da. 

JjQ*x(0,2r) JQ* 

This concludes the proof as 

n „ 

l|V«(.,0(.))||| 2(Q ) < / w^Of^dx. 

i=l ^ n_1 

□ 

From now on we follow the stopping time argument from [T5j, in particular our 
Lemma H~4l is an analogue of [T4"l Lemma 3.8]. For any continuous function v : R™ — > R n 
and v G R we define 

h u , a (y)(x) = sup{t > 0; sup |v| > u}. 

r a (x,t) 

Here T a (x,t) is a cone with vertex at (x,t) (recall that the boundary point is (x,0)). 
Hence 

r a (x,t) = (o,t) + r a (x,o), 

is the non-tangential cone T a (x, 0) shifted in the direction (0,t). 

Lemma 4.5. Ifv is such that h Uia (y) < oo then h Uja (v) is Lipschitz with constant 1/a. 

Proof. See, for example [21 Lemma 3.13]. □ 

We also have an analogue of [T4"l 3.14]. 

Lemma 4.6. Let u be a solution to Lu = divAVu = 0, where L is an elliptic dif- 
ferential operator defined on R™ with bounded coefficients which are such that A3.1\) is 
the density of a Carleson measure on all Carleson boxes of size at most ro with norm 
e. Set v = Vn and let (Qj)j be a Whitney decomposition of {x; iVj ](v)(a;) > ^/24}. 
Given a > 0, let E ] u be the intersection of the cube Qj with 

{x; N[ a / 12 ](v)(x) > v and EN[ a ](v)(x) + S[ a] (v)(x) < pv). 

Then there exist a sufficiently small choice of p, independent of Qj so that for each 
x G El there is a cube R with x G 6R and R C Q* for which 

|v0z,/l„, tt /12(v)(2))| > v/2 

for all z G R. 

Proof. Let x G El . By definition /^^ a /i2(v)(x) > and so there exists a Y on 
dT a/12 (x, h v>a / 12 (v)(x)) such that \v(Y)\ = v (here Y = (y,y n )) and h u>a / 12 {v){y) = y n - 
Let r = y n and 

K = r a/ i 2 (x,0) n {Z; \z n - y n \ < r /6}. 
Since Qj is a Whitney cube, r < (1 + A\/n — l)£(Qj)/a, and we also have 

3K C r o (z, 0) and 6M{3K, 9R+) > r /2. 
Hence again by [HI Thm 8.17] we have that 



BOUNDARY VALUE PROBLEMS FOR OPERATORS SATISFYING CARLESON CONDITION 25 



osc K (v) < C(r - /2 ||v - c\\ LH2K) + rl^UVAM^), 
for any constant c and q > n. By (13. ip \{V A)w\(Z) < CrQ l eN[ a ](v)(x) for Z e 2K, so 

r l - n,q \\(yA)v\\ Lq{2K) < CeN [a] (v)(x) 
and so using Poincare's inequality 

\v(Z)-v(Y)\ < osc^(v) <C(r;- n/2 ||Vv|| i2{3i ,)+£iV [a ](v)(x)) 

< C(S [a] (v){x) + eN [a] (v)(x)) < Cpv, 

for any Z G K. Thus we may choose p sufficiently small so that |v(Z) — v(Y)| < vj2. 
Then clearly \v(z, h Uta /i 2 (v)(z))\ > is/2 for \z — y\ < ar /72. □ 

Finally, the good-A inequality can be converted using standard methods to a global 
inequality between N and S. We have the following: 

Lemma 4.7. Let u be a solution to Lu = divAVu = 0, where L is an elliptic differ- 
ential operator defined on Q to with bounded coefficients such that 

(4.17) sup{5(X)|Va 4J (F)| 2 : Y e B S{X)/2 (X)} 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(r ). 

Then there exists e > depending only on the geometry of the domain Q, the ellip- 
ticity constant A, dimension n and p such that if C(r ) < e then 

(4.18) f N 2 r/2 (Vu)dx<K [ S*(Vu)dx+ [[ \Vu\ 2 dX. 
Jan Jan JJa\a r/2 

Here K only depends on the geometry of the domain Q, elliptic constant A, p and 
dimension n. Here N^, Sh are truncated versions of non-tangential maximal function 
and square function, respectively. 

Remark. The term ff n \ n / |Vw| 2 dX is necessary if Q is a bounded domain. Consider 

for example L = A on O C 1". Let u be a harmonic function in Vt. Then for any 
vector c we have that S(Vu) = S(V(u + c • x)) but clearly N(Vu) ^ N(V{u + c • x)). 
This term is not necessary if the domain is unbounded and we consider untruncated 
versions of the non-tangential maximal function and the square function. 

Proof. We only highlight the major points of the proof as the basic idea is the same as 
in [T4]. Applying standard good-A inequality methods and Lemma 14.6} we see that 

|{z;JV [o/13] (Vu)(a;) > A}| 

< I {x; N [a/12] (Vu) (x) > X, eN [a] (Vu) (x) + S [a] (Vu) (x) < p\} | 
+ \{x;eN [a] (Vu)(x) + S [a] (Vu)(x) > pX}\ 

< \{x; M(Vu(; h u>a/12 (Vu)(-))(x) > X/2}\ 
+\{x;eN [a] {Vu)(x) + S [a] (Vu)(x) > pX}\ 
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Multiplying this inequality by A and integrating in A, we obtain the inequality 

\\N [a /i2](v)\\ 2 L2 {Qo) 

< c (||Af( v (-,^a/i 2 ( v )(-)))lli 2 (Q„) + II^mC^III,^) + e||Ar M ( v )||| 2(Qo) ) 

(4.19) < C (||v(-,^ a/12 (v)(-))||i 2(Q) + ||S [a] (v)||i w + e\\N [a] (v)\\ 2 L2m ) . 

Here M is the Hardy-Littlewood maximal function, Qq is a cube, Q = 2Qq and Q* = 5Q 
is such that Q* is contained in a single coordinate patch Ug. 

Now we can apply Lemma 14.41 to the first term on the right-hand side of (14.191) and 
after summing over an appropriate collection of cubes Qo obtain 

< C (\\S [a] (Vu)\\ 

L 2 {dn) 

(4.20) + \\N [a] (Vu)\\ L 2 im) \\S [a] (Vu)\\ L 2 (m) + if \Vu\ 2 dX), 

JJn\n r/2 

Note that Lemma 14.41 requires the Lipschitz function <fi has a small Lipschitz norm. 
Since we are using function /i„, a /i2(v) in place of 0, if we choose a > large enough 
by Lemma 14.51 the Lipschitz norm will be small. Standard techniques also tell us that 
\ \N[ a ]( Vu)\\ L 2 {dn) < ||iV[ a /i2](Vu)||i,2( 9n ), S o if s is chosen small enough in (14.20}) then 
(|4TT8|) will follow. □ 

5. The (R) 2 Regularity Problem 

Theorem 5.1. Let Q C M be a Lipschitz domain with Lipschitz norm L on a smooth 
Riemannian manifold M and Lu = div(AVu) be an elliptic differential operator defined 
on fl with ellipticity constant A and coefficients which are such that 

(5.1) snp{6(X)\V aiJ (Y)\ 2 : Y e B 5{X)/2 {X)} 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C{tq). Then there exists e = e(A, n) > such that if max{L, C(r )} < e then the 
regularity problem 

Lu = 0, in f2, 
u = f, on dfl, 
N(Vu) G L 2 (dn), 

is solvable for all f with \\ Vt/Hl 2 ^) < oo. Moreover, there exists a constant C = 
C(A,n,a) > such that 

(5.2) HA^WllU^a) < C||V T /|U 2(9n ). 

Proof. For any / in the Besov space B 2,2 2 (dQ) the exists a unique H 2 (Q) solution by 
the Lax-Milgram theorem. Observe that / £ H 2 (dQ) C B 2 ', 2 (dQ) so it only remains 
to establish the estimate (15. 2p . 

Consider e > and take C(tq) < e. To keep matters simple let us first consider the 
case when dQ is smooth. In this case Lemma 13.21 applies directly. If follows that for 
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all small r 



(ay 



(5.3) f S 2 r JVu)da<K [ |V T u\ 2 da + C \e) [ N 2 (Vu)da + — ||Vu||£a 
Jan Jan Jan r 

Here C(e) — > if e — > 0. We now choose e small enough such that Lemma [4.71 holds. 

It follows that 

(5.4) / N 2 /A {Vu)da<K [ \V T u\ 2 da + C{e)K 2 [ N 2 {Vu)da 
Jan Jan Jan 

+ K(r)\\Vu\\ 2 LHn) . 

Here K appearing in ( 15 .4p is the same constant as in the estimate ( I4.18p . 
We also observe that we have a pointwise estimate 

(5.5) N 2 {Vu){X)<N 2 u{Vu){X) + C{r)(( \Vu{Y)\ 2 dY 



for all X G dfl. Seeing this is not hard as we are estimating | Vm| away from the bound- 
ary. Hence by the Carleson condition we have |Vy1| < C(r) there. Rest is a standard 
bootstrap argument using the equation v = V« satisfies, i.e., Lv = div((VA)v) even- 
tually yielding pointwise bound on |Vu| for {X G dQ; dist(X, dQ) G [r/4, r]}. Finally, 
combining (15.41) and (15. 5p we obtain 

(5.6) f N 2 (Vu)da<K ( \V T u\ 2 da + C(e)K 2 [ N 2 (Vu)da 
Jan Jan Jan 

+ k(r)\\Vu\\ 2 L2{n) . 

We now can make our final choice of e. We choose it sufficiently small such that the 
constant in (15. 6p C(e)K 2 < 1/2 which yields 

(5.7) f N 2 (Vu)da<2K [ \V T u\ 2 da +2K{r)\\ Vu||2 
Jan Jan 

From this the desired estimate follows since the term || Vm||| 2( -q), i-e., an Hf(Q) estimate 

of the solution u, follows from the Lax-Milgram. 

Now turn to the more general case, when Q has a Lipschitz boundary with sufficiently 
small Lipschitz constant L. This case also includes the C 1 boundary as in such case L 
can be taken arbitrary small. 

The crucial point is that the proofs of Lemmas I3.2H4.7I in the smooth case are based 
on local estimates near boundary dQ. This means we can reduce the matter to a 
situation where we have an open set U in M" and a Lipschitz function with Lipschitz 
constant L such that in U the set Q looks like {(x, t) G M n ;t > (f>(x)}. 

Now the map (14.101) is a bijection between the sets M™ and {(x,t) G M. n ;t > <f)(x)}. 
In fact if c > £ then the map $ is a local bijection,, where £ = \\V(/)\\bmo- Hence 
by pulling back everything (metric, coefficients) using $ we are left with proving local 
estimates on a subset of 1R™. However, this is exactly what we did above. We only 
have to be careful about how much the Carleson constant of the coefficients changes 
when we move from the set {(x, t) G M n ; t > <f)(x)} to 1R™ . A computation gives us that 
if the original constant was C, the new constant on will be C + C(£) where C(£) 
is an increasing function in £ such that lim^ + C(£) = 0. From this the claim follows, 
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as this implies that C + C(£) will be small as long as both C and i are small enough. 
So we get solvability on domains with small Lipschitz constant, as well as on domains 
whose boundaries are given locally by functions with gradient in VMO. □ 

Finally, we replace the gradient Carleson condition ( 15.1 \ by a weaker condition for 
oscillation of the coefficients ( 15. 8p . What this means is that the gradient Vw no longer 
has well-defined pointwise non-tangential maximal function N. Instead a weaker ver- 
sion N defined by (12.11) must be used. 

Theorem 5.2. Under the same assumptions as in Theorem I5.il the (i?) 2 regularity 
problem for the operator Lu = div(AVu) is solvable under a weaker Carleson condition 
requiring that 

(5.8) ^W" 1 {osc B{ x,s(x)/2)aij) 2 

is the density of a Carleson measure on all Carleson boxes of size at most r with norm 
C(r ) < e. 

Moreover, there exists a constant C = C (A, n, a) > such that 

(5.9) \\N(Vu)\\ L2{m) <C\\V T f\\ L 2 m . 

Proof. The proof uses same idea as [5j Corollary 2.3]. For this reason we skip the non- 
essential details. The procedure outlined in [5] implies that for a matrix A satisfying 
(15. 8p with ellipticity constant A one can find (by mollifying coefficients of A) a new 
matrix A with same ellipticity constant A such that A satisfies (15 .ip and 

(5.10) supWXy'KA- A)(Y)\ 2 ;Y e B(X,5(X)/2)} 

is the density of a Carleson measure. Moreover, if the Carleson norm for matrix 
A is small (on balls of radius < r ), so are the Carleson norms of (15.11) for A and 
(I5.10p . Hence by Theorem 15.11 the (i?) 2 regularity problem is solvable for the operator 
Lu = div(AVu). 

The solvability of the regularity problem for perturbed operators satisfying ( I5.10P 
has been studied in [13]. It follows by [T3| Theorem 2.1] that the LP regularity problem 
for the operator L is solvable for some p > 1. The p for which the solvability of the 
regularity problem is assessed is the p such that the LP , p' = p/(p — 1), Dirichlet 
problem for the adjoint operator L* is solvable. Actually, the results in [TH] are stated 
for symmetric operators, however careful study of the proof of [T3"| Theorem 2.1] reveals 
that what is really needed is to replace L by its adjoint when the LP Dirichlet problem 
is considered. 

However by [5j Theorem 2.2] the L 2 Dirichlet problem for L* is solvable provided 
the Carleson norm of (15. 8p (and hence (15. ip for A) is sufficiently small. Hence we have 
solvability of the regularity problem (i?) 2 by [131 Theorem 2.1, Remark 2.3]. □ 

6. The Square Function Revisited 

In this section we revisit bounds for the square function of Vw from the perspective 
of the Neumann problem. As in Section 2 we shall assume that Q is a smooth domain 
on a smooth compact Riemannian manifold M. We also continue to use notation we 
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introduced there. Recall that Q to denotes the collar neighborhood of the boundary 
dQ x (0,t ). 

On Q to we have a well-defined co-normal derivative of u with respect to the operator 
L; in the metric da £*D dt this is just 

i=l 

where (ay) are coefficients of the matrix A in local coordinates near the boundary. 

We have the following key lemma bounding the non-tangential maximal function of 
Vw by the square function of H. 

Lemma 6.1. Let p > 2 and u be a solution to Lu = divAVu = 0, where L is an 
elliptic differential operator defined on a smooth domain Q with bounded coefficients 
such that M-lty is the density of a Carleson measure on all Carleson boxes of size at 
most r with norm C(r ). Then there exists e > such that if C(r ) < e then for some 
K = K{Q,A,n,e,p) > 

(6.1) f N p {Vu)dx<K ff \V T u\ p - 2 \VH\ 2 5(X)dX + C(r) ff \Vu\ p dX. 
Jan JJn 2r JJn\n r/2 

Proof. We mainly work in the collar neighborhood Q to defined above. We choose r < 
to/5. Using the results we have on the solvability of the regularity problem we know 
that for sufficiently small e > we have: 

(6.2) f N p {Vu)dx<K f \V T u\ p dx. 

Jan Jan 

Since dQ, is a smooth compact manifold, there is a finite collection of balls Qx, Q%, . . . , 

in M n_1 of diameter comparable to r and smooth diffeomorphisms ipe : 5Qp —> dfl such 

that [j e ipi(9/8Qt) covers dfl. Here rQ denotes the concentric enlargement of Q by a 

factor of r. Let us also find smooth partition of unity <p£ such 

<fii = 1 on dfl, (fie = 1 on Q t and supp (fii C 9/8Qe. 

Let us fix i and work on one ball Q = Qe and £i = (fig. We may assume that 
Wi\ — C/r. Choose another function £2 : [0, 00) — > R such that ^{t) = 1 on [0,r], 
1^2 1 — 5/r and support contained in [0,2r]. Now define 

(6.3) £(X) = £0M) = 6 

We work on estimating righthand side of (16.21) in local coordinates on 5Q x (0,5r). 
Denote by Vk = d^u for k = 1, 2, . . . , n. For each k < n — 1 we have 

f \v k (x,0)ni(x)dx = - ff d n (\v k \ P 0(X)dX 

(6.4) = -p ff \v k \ p - 2 v k (d n v k )£dX- ff \v k \ P W 2 dX = I + II. 

The second term on the right-hand side of (16.41) is controlled by JJ K |Vu| p where 
K = {X = (x,t);x G 5Q and r/2 < t < 5r}. We deal with the first term. Since 
d n v k = d k v n we have 
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P // \v k \ p - 2 v k (d k v n )£dX 



(6.5) = -p 



[[ wr 2 v k d k (**v\ idx+p^ff wr\ k d k (**v\ tdx. 

JjR 1 } V^nn / . JJr" \ a nn J 



i<n + 

The second term of (16.51) can be further written as 



pY,IJ \v k r 2 v k d k (^v)tdx 
(6.6) =pY / [[ i^rv^N^i+Ef dMn—tdx. 

We introduce (d n t) into both the terms of ( 16.61) and integrate by parts. This gives 



P 



[[ d n (\v k r 2 v kVl d k AtdX + [[ d n (diUvtf)^) tdX 

J J R™ V \ a nn J J JjWL V a nn J 



-E 

i<n 

-E 

i<n 

(6.7) + // di(\v k \*)d n (^)&dX 



P // d n (\v k r 2 v k ) Vi d k [^L)^tdX + p \v k r 2 v k d n ( Vi )d k (^)itdX 



E 



E 



p \v k r 2 v k v l d k (^)w 2 tdx+ dl (\v k n—ti&dx 



P II \v k r 2 v kVl d n d k (^\£tdx+ dndi (\ Vk \p)^tdx 

&nn J JjR" a nn 



The last two terms we integrate by parts one more time as we switch the order of 
derivatives. This gives 



R™ V a nn 



E p II d k {\v k r 2 v k v£)d n (^)tdx+ ll ^(^)«9 n (kntdx 

i<n 

= E 

i<n 

(6.8)+ // d n (\v k \ p )d i [ — )CtdX 



P // d k (\v k r 2 v k )v t d n ^ itdX + p // l^r^C^i;*)^ (^)itdX 



Kn 



P II \ v k\ p 2 v k Vid n 



(d k ^ 2 tdX+ I/ d n (\v k \ p ) — {di^ 2 tdX 
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The first three terms on the righthand side of both (16.71) and (16. 8p can be bounded 
from above by 



|Vtt| p_1 |V 2 tt||VA|^X 

2Qx[0,2r] 

1/2 / rr x 1/2 



(6.9)C ff 

JJ2i 

< ( ff \Vu\ p - 2 \\7 2 u\Hdx) ' ( ff \Wu\ p \WA\ 2 tdx) 

\Jj2Qx[0,2r] J \Jj2Qx[0,2r] J 

f N(Vu) p ~ 2 ff \V 2 u{X)\ 2 t 2 ~ n dX dx) ' ( ff \Wu\ p \VA\ 2 tdx] ' 

J2Q JJr(x) J \JJ2Qx[0,2r] J 

f N{Vuy- 2 S 2 {Vu)dxX'\\\N{Vu)\\% Q ^ 



(6.10) C I \Vu\ p \VA\±dX 

'2Qx[r,2r] 



< ( ff \Vu\ p ± 2 dx)' ( ff \\7u\ p \X7A\ 2 tdxY 

\Jj2Qx[r,2r] J \Jj2Qx[0,2r] 

^ N(Vu) p (x) dx) e\\N(Vu)\\ P l\ 2Q) = e\\N(Vu 



2Qx[r,2r] J \JJ2Qx[0,2r] 

1/2 

V 

Lp(2Q)- 

'2Q / 

The fifth term on righthand side of (16.7p can be estimated by 
C ff \Vu\ p ^\V 2 u\^ r dX 

JJ2Qx[r,2r] 

/ rr \ P/(P-1) / rr \ 1/p 

(6.11) < // \Vu\ p dX) // \V 2 u\ p dX 

\Jj2Qx[r,2r] J \Jj2Qx[0,2r] 

< C{r) JJ \Vu\ p dX. 

To get the last line we used some standard elliptic estimates away from the boundary 
(for example, it is sufficient to generalise Caccioppoli's inequality to inhomogeneous 
equations via the proof in [10], p. 2]). By the Carleson condition we have | VA\ < C(r) 
there. Rest is a standard bootstrap argument using the equation v = Vu satisfies, i.e., 
Lv = div((VA)v) eventually yielding L p bounds on Vv in K. 

We denote the co-normal derivative of u by H = a n idiU = ^ a n (Vi and write the 
first term of (16. 5p as 
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-V ff \v k \ p - 2 v k d k (—) £dX = -p ff \v k \ p - 2 v k d k (—) ad n t) dX 
= P ff d n (\v k r 2 v k )d k (—)^tdX + pff \v k \ p - 2 v k d k (—)W 2 tdX 

JjR™ \ a nn J JjR™ \ a nnJ 

(6.12) +p ff \v k \v- 2 v k d n d k (—)&dX, 

JJM.™ V a nn/ 

where the last term further yields: 

(6.13) -p ff d k (\v k r 2 v k )d n (—) ^tdX-p ff \v k r 2 v k d n (—) (d k ^ 2 tdX. 

JjR" \ a nnJ JJm. 7 } \ a nnJ 

If the derivative in the first two terms on the righthand side of (16. 12)) and ( 16. 13ft falls 
on the coefficients of the matrix A we obtain terms we have already bounded above 
(see (I6.9P and ( 16 . 1 0[) ) . If the derivative falls on H the first term on the righthand side 
of both (16U2|) and (T6TT3]) is bounded by 

C ff \v k \ p - 2 \Vv k \\VH\&dX 



V2 / \ 1/2 

< C[ II \v k \ p - 2 \Vv k \ 2 &dX) // \v k \ p ~ 2 \VH\ 2 itdX 



< C 



1/2 > \ V2 



f N p - 2 (v k )(x) f \Vv k (X)\H 2 - n dX dx] ( ff \V T u\ p - 2 \VH\ 2 &dX 

J2Q JT(x) J VJJm^ 

C^N p - 2 (v k )(x)S 2 (v k )(x)d x y (jj^ \V T u\ p ~ 2 \VH\ 2 itdX 



'2Q 

If/ 2 " 1 II .<?/'„. Ml I II |V7„,,|P-2| 



= CHJV^II^jII^Hlpw) I jj^ \V T ur 2 \VH\ 2 itdX 

If the derivative falls on H in the second term of (16.121) . we get terms of the same 
form as doUOj) and ( IBTTTjl . 

It follows that for all k < n — 1 we have 



(6.14) / v k (x,0) p ^(x)dx 

JR"- 1 

< ^11^(^)11^)11^(^)11^(20) ( £ n \V T ur 2 \VH\ 2 itdX\ 

+e||5(V«)||xP( a0 )||MV«)||^J l g ) +s\\N(Vu)\\ p LP{2Q) + C(r) JJjVu^dX 
+E. 
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Here E denotes remainder terms; these are the last two terms of (16. 8p and the last term 
of (I6.13P when the derivative falls on H. We now sum (16.141) over all k < n — 1 and also 
sum over all coordinate patches Qg. We notice that the error terms E with complete 
cancel out as \\^ t {dk<pe) = where (4>e)e is the partition of unity we considered above. 
This yields a global estimate 

\V T u\ p dx < ^11^)11^115(^)11^) ( [[ \V T ur 2 \VH\ 2 5(X)dx) ' 
an \JJn 2r J 

+e\\S{Vu)\\ Lvm \\N{Vu)f-^ +e\\N(Vu) \\ p LP{da) + C(r) [[ |V<dX 

JJn\n r/2 

From this, by (16. 2p and using Lemma I4.3[ we get that for all sufficiently small e > 
the desired estimate (16. ip holds. □ 

Lemma 6.2. Let p > 2 be an integer, k be an integer such that < k < p — 2 and u be 
a solution to Lu = div AVu = 0, where L is an elliptic differential operator defined on a 
smooth domain Q with bounded coefficients which are such that jjj4.1 r /\ ) is the density of a 
Carleson measure on all Carleson boxes of size at mostro with norm C(ro). Then there 
exists e > such that if C(ro) < e then for some constant K = K(Q, A, n, e, m,k) > 



(6.15) [[ \V T u\ p - k - 2 \H\ k \VH\ 2 5{X)dX 
JJn r 



< {p-k-2)K \V T u\ p ~ k - 3 \H\ k+1 \VH\ 2 5(X)dX + C(r) \Vu\ p dX 
JJn 2r JJn\n r 

+K [ \H\ p dx. 
Jan 

Proof. We will establish (I6.15P by induction on k. If k = by Lemma [6. II we have: 



(6.16) / N p (Vu)dx<K \V T u\ p - 2 \VH\ 2 5(X)dX + C(r) \Vu\ p dX. 
Jan JJn r JJn\n r/2 

For k > we use (16.11) and the induction assumption (I6.15P for all indices 0, 1, . . . , fc — 
1. This gives 

f N p (yu)dx < K [[ \V T u\ p - k - 2 \H\ k \VH\ 2 5{X)dX + C{r) [[ \Vu\ p dX 
Jan JJn r JJn\n r/2 



(6.17) +K / \H\ p da. 

Jan 

Here K = K(k) and (I6.17P holds for all sufficiently small s > 0. From this, the 
inequality 
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/ N p (Vu)dx + K ff \V T u\ p - k - 2 \H\ k \VH\ 2 S(X) dX 

(6.18) < 2K ff \V T u\ p - k - 2 \H\ k \VH\ 2 5{X)dX 

+2C(r) ff \Vu\ p dX + 2K f \H\ P da 
JJn\n r/2 Jan 

holds when k = without any further assumptions, and when k > under the induc- 
tion hypotheses (16.151) for indices 0, 1, . . . , k — 1. Let us choose a cutoff function £ as 
in (16.31) . To control 



ff \V T u\ p - k - 2 \H\ k \VH\ 2 5{X)dX 

J J fir 



it suffices to control 



ff \V T u\ p - k - 2 \H\ k b ij (d i H)(d j H)£tdX = I 

for some matrix B satisfying the ellipticity condition to be specified later. 
We integrate this by parts. This gives 

7 = ff \V T u\ p - k - 2 \H\ k Hd i (b ij d j H)^tdX 



TXT ff \V T ur k ~ 2 \H\ k Hb nj {d 3 H)idX 
& + t JJw\ 



(6.19) ff iVTur^lHl'Hb^HMOtdX 

- P ~ k ~ 2 ff \V T u\^ k ^(V 7 ^-d i (V T u))\H\ k Hb nj (d j H)&dX. 

The second term only appears in (I6.19P if z = n as the function i obviously only depends 
on the variable x n = t. We first deal with the third term of (16.191) when i = n. As 
| £2 1 < 2/r and £ 2 = on [0, r] we have that this term is bounded by 

(6.20) /" / |V T u| p - fe - 2 |iJ| fc+1 |ViJ|^X < £ f N p r (Vu)dx + C{r) ff \Vu\ p dX, 

JQ Jr J2Q JJn\Q r 

since this term is of same type as (I6.10P and (16. lip it can be estimated as before. 
Now for the terms with % < n in the third term of ( I6.19P we observe that they will 
cancel when we sum over the index I in the partition of unity introduced via the cutoff 
function £ from (I6.3p . We choose the matrix B so that b nn = 1. Then the second term 
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of (I6.19P if j = n looks like 

1 



(k + l)(k + 2) 
1 

'(fc + l)(fc + 2) 



\V T u\ p ~ k ~ 2 {d n \H\ k+2 )^dX 
d n {\V T u\ p - k - 2 \H\ k+2 i)dX 



+ _L // \V TU \ p - k - 2 \H\ k+2 i'dX. 



{h + l){k + 2) 

Here the last term again can be estimated by a solid integral C(r) JJ n ^ n \Vu\ p dX in 
the interior of the domain. The first term is equal to a boundary integral 

(jfc + 1)^ + 2) L WTur k - 2 \H\ k + 2 dX < v\\^Tu\\ p LP{dn) + C(rj)\\H\\l P{m) , 

for rj > arbitrary small. Note that 

^l|VTn|r LP(M) <r/||iV(Vn)r LP(9n) . 

We choose rj > so small that we can hide the term v\\N(Vu)\\ p LP ( dn - ) on lefthand side 
of (J6H]). 

It remains to deal with the second term of (I6.2ip . We differentiate and change the 
order of derivatives d n and Vt- 

— 2 ff iv-7 „.io-fc-4/T-7 „, vy a „.\\TT\k+2i 



(6 ' 22) (* +!)(* + 2) JL ™'-*-\Vru • Wm^i dX . 

We reintroduce the co-normal derivative H as d n u = — — —Vj. We also insert 

Q>nn — J < ~~''- (Inn J 

a term (d n t) = 1 into both integrals. Then we integrate by parts again in the d n 
derivative. Whenever exactly one derivative falls on the coefficients (either a nn or ^ 
those terms are bounded by 



Tin ' 



(6.23) ff 

JJ2C 



|V^||Vu| p_1 |V 2 M|tdA' 



'2Qx[0,2r] 

which is the term of type ( 16. 9p and has therefore a bound of type £||iS'(Vu)||lp(2Q)||A 7 "(Vw) 
with e depending on the Carleson norm. For sufficiently small e, thanks to Lemma fl~3l 
this can be hidden on the lefthand side of ( 16. 18ft . 

If both d n and Vt derivative fall on coefficients, there are two possibilities. The 
first possibility is that they fall on the same coefficient and so then we do a further 
integration by parts in Vt moving this derivative on other terms. This again will yield 
term of type (16. 23[) . The second possibility is that they fall on separate coefficients 
and so take the form ( I4.16p . which can be estimated appropriately with the help of 
Lemma 14.31 We obtain another error term when d n falls on £, however in that case 
we get a term of type ( I6.20p we handled before. Let us deal with the term when both 



liP(2Q)' 
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derivatives fall on H. In that case we have 

(6.24) - ./~f.~ 2 . ff —\V T u\ p - k -\V T u-V T d n H)\H\ k+2 ZtdX. 



(k + l)(k + 2) JJ n2r a n .. 

We move the Vt derivative off d n H . We can get a term of type (16.231) and two terms 
that can be dominated by 

(6.25) Op-k-2) fl \V T u\ p - k - 4 \V(V T u)\\VH\\H\ k+2 tdX 

+ Op-k-2) ff \V T u\ p ~ k ^\VH\ 2 \H\ k+ HdX. 

Also, when Vr lands on £ we get error terms which will cancel when we sum over 
coordinate patches. Observe also that the last term of ( 16. 19f) can be controlled by 




(6.26) C(p-k-2) // \V T u\ p - k ~ 3 \V(y T u)\\VH\\H\ k+1 tdX 

We now deal with the terms arising from — V . , —v^. Here we write 

= E v -(? i )^+E^( v -")- 

The contribution of the first term here, when substituted in (I6.22p . can be dealt with 
by again introducing the factor d n t and integrating by parts. When d n lands on 
^Tidnj/O'nn), w e can move the tangential derivates off by again integrating by parts. 
All this yields terms of the form (I4.16P (with exponent p instead of 2) and (I6.23p . which 
can be controlled appropriately. Substituting the second term in (I6.22p yields 

J<~.Tl 

Moving d n across using integration by parts and if necessary moving dj we obtain terms 
either bounded by (14.161) (with exponent p instead of 2), (16.231) . (I6.25P or (16. 26[) . Thus 
the analysis of the second term of (16.191) for j = n reduces to controlling (16.251) and 
(I6.26p . a task which we will postpone for now. When j < n in the second term of 
(I6.19P we again introduce (d n t). This gives 



II \V T u\ p ^ k ~ 2 b nj dj(\H\ k+2 )^(d n t) dX 



We integrate by parts. When d n falls on \W Tu\ p ~ k ~ 2 we can dominate such a term by 
(I6.26p . when d n falls on b n j we obtain a terms of type (14.161) (with exponent p instead 
of 2) and (I6.23P and, provided we choose matrix B so that coefficients of B also satisfy 
the Carleson condition. If d n hits £ we get terms which can be bounded by (I6.10p and 
(16. lip . Finally the remaining term is 

ff \v T u\ p - k ~ 2 b n3 d 3 d n {\H\ k+2 )itdx. 

We integrate by parts again in dj giving us terms of type (I4.16P (with exponent p instead 
of 2), (I6.23p . (I6.25P and (16.261) . The only remaining terms we have not yet bounded 
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are the first term of (16.191) . (16.251) and (16.261) . The second term of (I6.25P is already 
of desired form (see righthand side of (I6.15P ). By the Cauchy-Schwarz inequality, the 
first term of (16.251) can be bounded by 

1/2 



Cip-k-2) (^jj \V T u\ p - k - 3 \VH\ 2 \H\ k+1 tdX^J x 



2r 

1/2 

\V T u\ p - k - 5 \V(V T u)\ 2 \H\ k+3 t dX 

1/2 



< C(p - k - 2) 
The last line can be further bounded by 



£ \V T ur k - 3 \VH\ 2 \H\ k+ hdX^j ||iV(V M )|r L / p 2 ( ^ ) ||5(V M )|U P( an ) . 



v\\N(Vu)\\% m + C( V )(p-k-2)* If \V T ur k ~ 3 \VH\ 2 \H\ k +HdX, 

for 7] > arbitrary small. Hence as before we can hide ^HA^Vm)!^^ on the lefthand 
side of (16.181) . Term (I6.26P can be dealt with in a very similar fashion. We summarize 
what we have so far. By (16 . 1 8[) and all estimates above we have 



a 



mi 



N p {Vu)da + [[ \V T u\ p ~ k ~ 2 \H\ k \VH\ 2 5{X)dX 
JJn r 



(6.28) < K(p-k-2) // \V T u\ p ~ k - 3 \H\ k+1 \VH\ 2 8(X)dX 



2C(r) \l \Vu\ p dX + K I \H\ p da 
JJn\n r/2 Jan 

— — — If \V T u\ p - k ~ 2 \H\ k H{LH)&dX. 
k + 1 JJnor 



Here LH = div(B'VH) and a > 0. The precise value of a depends on choice of rj > 
above and e > 0. Clearly, (I6.28P is the desired estimate (I6.15P modulo the last extra 
term we shall consider now. 

As above we use the summation convention, we only write the sum explicitly when- 
ever we do not sum over all indices. For LH we have 

LH = ditbijdjH) = ^ditbijdjtcinkdku)) + di(b in d n (a nk d k u)) . 

j<n 

Since Lu = we have that d n (a n kdku) = — J2j< n ^j( a j^ku). Hence 

LH = di(bijdjH) = S ^2 / [d i {b ij d j {a nk d k u)) - di(b in d j (a jk d k u))]. 

j<n 

We also swap the role of % and k in the second term. From this 

(6.29) LH = diibijdjH) = ^2[d i (b ij d j (a nk d k u)) - ^(^^(a^n))]. 

j<n 
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We choose b^ = aji/a nn . Notice that this guarantees that b nn = 1 and that terms in 
( I6.29P where three derivatives fall on u vanish as these are the terms: 

(6.30) ^[bijan^didjdku) - 6 fcn aj i (9 i 9,-9 fc u)] = ^a^a^a^ - a^a^did^u = 0. 

j<n j<n 

We now place (I6.29P into last term of (I6.28p . Given (I6.30p some of the remaining 
terms are 



(6.31) ff Yihjidid^iduu) -bknid^a^id^WVTU^-^H^H^tdX 
JJn 2r j<n 



and the rest can be bounded by 



(6.32) // \Vu\ p - L [\Vu\\VA\\VB\ + \V 2 u\\VA\\B\ + \V 2 u\\VB\\A\]tdX. 

The terms in (16.311) have two derivatives on coefficients however one is dj and j < n. 
We therefore integrate by parts in dj. This yields additional terms, but all are of the 
form (16.321) . However, by an estimate similar to (16.231) we get that all the terms of 
(I6.32p are smaller than e J dn N p r (Vu) da, with e depending on the Carleson norm of 
the coefficients of our operator. Hence for sufficiently small e this term can be hidden 
in (16.281) within the term a J dn N p (Vu) dx. This yields the desired estimate (I6.15p . □ 



7. The (N) p Neumann Problem 

Theorem 7.1. Let p > 2 be an integer and let Q C M be a Lipschitz domain with 
Lipschitz norm L on a smooth Riemannian manifold M and Lu = div(AVu) be an 
elliptic differential operator defined on Q with ellipticity constant A and coefficients 
such that 

(7.1) sup^pOlVa^Om 2 : Y £ B 5[X)/2 {X)} 

is the density of a Carleson measure on all Carleson boxes of size at most ro with norm 
C(ro). Then there exists e = e(A,n,p) > such that if max{L, C(r )} < s then the 
Neumann problem 

Lu = 0, in Q, 

AVu ■ v = f, on <9f2, 
N(Vu) £ L p (dn), 

is solvable for all f in L p (dQ) < oo such that J dn fda = 0. Moreover, there exists a 
constant C = C(A,n,a,p) > such that 

(7.2) \\N(Vu)\\ LP{d ^<C\\f\\ LP{d ^. 



Proof. For any / in the Besov space -B^ 2 (<9f2) such that J dn fda = the exists a 
unique (up to a constant) Hf(Q) solution by the Lax-Milgram theorem. Observe that 
our / £ L p (dVl) C B^y 2 (d£l) (p > 2) so it only remains to establish the estimate (17.21) . 
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P 

lp (n\n r )- 



Consider e > and take C(r ) < s. To keep matters simple let us first consider the 
case when dQ is smooth. In this case Lemmas 16.11 and 16.21 apply directly. If follows 
that for all small r and e > 

(7.3) / N p (Vu)da<K [ \AVu ■ v\ p da + C(r) || Vu\ 
Jan Jan 

Here we are using Lemma [6.21 for k = 0,1,2, ... ,p — 2 while observing that for the 

integer k = p — 2, the first term on the righthand side of (16.151) is zero. As AVu ■ v = f 

we have for non-tangential maximal function 

(7.4) f N*>(Vu)da<K [ \f\ p da + C(r)\\Vu\\l 
Jan Jan 

We also observe that we have a pointwise estimates on Vu(X) for all X away from 

the boundary. There, by the Carleson condition, we have |VA| < C(r). The rest 

is a standard bootstrap argument using the equation v = Vw satisfies, i.e., Lv = 

div((VA)v) eventually yielding pointwise bound on |Vw| for {X £ dQ; dist(X, dfl) > 

r}. 

This yields 

(7-5) HVti||^ (IAnp) < C(p)\\u\\ P Hm < C(p)\\f\\ P B ^ m - 

Finally, combining (j7.4p and (17.51) we obtain the desired estimate (17.21) . 

Now we turn to the more general case, when Q has a Lipschitz boundary with 
sufficiently small Lipschitz constant L. This case also includes the C l boundary as in 
this case L can be taken arbitrary small. 

The crucial point is that the proofs of Lemmas 14. 7\ 16.11 and 16.21 in the smooth case 
are based on local estimates near boundary dQ. This means we can reduce the matter 
to a local situation working on the subset of upper half-space via the map (I4.10p as we 
did above. From this the claim follows on domains with small Lipschitz constant, as 
well as on domains whose boundaries are given locally by functions with gradient in 
VMO. □ 
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